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Abstract 

We investigate an index theorem for a Bogoliubov-de Gennes Hamiltonian (BdGH) de- 
scribing a topological superconductor with Yang-Mills-Higgs couplings in arbitrary dimen- 
sions. We find that the index of the BdGH is determined solely by the asymptotic behavior 
of the Higgs fields and is independent of the gauge fields. It can be nonvanishing if the 
dimensionality of the order parameter space is equal to the spatial dimensions. In the pres- 
ence of point defects there appear localized zero energy states at the defects. Consistency of 
the index with the existence of zero energy bound states is examined explicitly in a vortex 
background in two dimensions and in a monopole background in three dimensions. 



1 Introduction 



Yang-Mills-Higgs systems admit topologically nontrivial field configurations such as vortices 
and magnetic monopoles [HIS]- They are considered to be responsible for various nonpertur- 
bative effects in particle physics. Recently topological insulators and superconductors have 
attracted much interest as a new phase of materials in condensed matter physics. Topological 
insulators are band insulators whose ground states are characterized by topological numbers 
[3]. The method of topologically classifying insulating ground states has been generalized 
to superconducting states described by the mean-field Bogoliubov-de Gennes Hamiltonian 
(BdGH). This provides a unified framework of the classification of non-interacting fermion 
systems with respect to time reversal and particle-hole symmetries [H El El E] • The topologi- 
cal classification was further extended to systems with topological defects [8] , which enables 
us to classify zero modes along line defects and zero energy bound states localized at point 
defects from the point of view of the topological universality class. 

For a superconductor with time- reversal symmetry (class BDI), there exists a unitary 
operator that anticommutes with the BdGH. It defines an extended chiral symmetry. Gener- 
ically, chiral symmetry ensures the topological stability of the zero energy states, since these 
states are controlled by the index theorem [9l [lOl [H] for the BdGH, implying the robustness 
under continuous deformations of the order parameters. Remarkably, the extended chiral 
symmetry can be defined in arbitrary dimensions. This is a sharp contrast to the usual chiral 
invariance in massless Dirac theories, which only exists in even dimensions. The extended 
chiral invariance also incorporates internal spin symmetry. By gauging the internal spin 
symmetry we can reformulate the system as a Yang-Mills-Higgs system, where the space 
varying gap parameters can be regarded as Higgs fields. 

In this paper we investigate the index theorem of the Dirac-type BdGH corresponding to 
the Yang-Mills-Higgs systems in arbitrary dimensions with special emphasis on the extended 
chiral invariance. Previously, zero energy bound states in two and three dimensions were 
investigated by Jackiw and Rebbi [12], and Jackiw and Rossi [13]. The index theoretical 
approaches of their models were explored by Callias [9] and Weinberg [10] . In particular, the 
Weinberg index theorem turned out to be quite useful for much more complicated systems 
such as non-Abelian vortices in a color superconductor [M]. We show that for more generic 
models with Yang-Mills-Higgs couplings in arbitrary dimensions the Weinberg index theorem 
is applicable, and the index of the BdGH can be nonvanishing, implying that there appear 
topological zero energy bound states localized at the point defect if the dimensions of order 
parameter space is equal to the spatial dimensions. We also see that the index is solely 
determined by the behaviors of the Higgs fields at the spatial infinities. This is expected 
in odd dimensions [9], since we have no topological invariant defined by the gauge field. 
In even dimensions, however, we have two topological invariants, one defined by the Higgs 
fields and one by the Yang-Mills fields. The expression of the topological invariants seems 
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to contain both of them. One might consider that the number of zero modes would be 
affected by including the Yang-Mills fields. We show that this is not the case. These two 
topological invariants are so combined that a unique index is obtained. This was noted for 
vortex background in two dimensions [10]. We show this happens in arbitrary dimensions. 

This paper is organized as follows. In Sect. [2] we introduce BdGH in arbitrary dimensions 
and their generalized chiral symmetry. We then define the generalized index for the BdGH. 
A computation of the topological index is given in Sect. [31 The chiral current at the spatial 
infinities is investigated in Sect. |H General properties of the chiral current and the relation 
with the topological index in even dimensions are given in Sect. |5l Zero modes for vortex and 
monopole configurations in two and three dimensions are investigated in Sect. El Finally, 
Sect. [71 is devoted to summary and discussions. 

2 Extended chiral symmetry of BdGH 

We begin with a d dimensional fermion system coupled to an 0{d) Higgs field (pa (a = 
1, • • ■ ,d). We assume that the system is described by the following Hamiltonian 

= -^7^'9,+0(x), (0(x) = r>,(x)) (2.1) 

where 7-' (j = 1, ■ ■ ■ , (i) and 7"^+" = r° form a set of 2d dimensional 7 matrices satisfying 
{7'^,7'^} = 5^'^ . We introduce the 2d dimensional chiral matrix by 72^+1 = (— z)"'7^ ■ ■ ■7^'^ = 
{—iY'-y^ ■ ■ ■7'^r-^ ■ ■ ■ V^. It anti-commutes with both 7-' and F". The system then possesses 
the chiral symmetry in the sense that 

{72d+i,i/o} = 0. (2.2) 

Note that the chiral symmetry can be defined in any dimensions. This is contrasted with 
the usual chiral symmetry which is only defined in even dimensions. The Hamiltonian 02. ip 
concerns itself with 2d dimensions, d spatial and d internal. The chiral invariance is known to 
be violated in the presence of the chemical potential. In the present work we are concerned 
with the chiral symmetric case, where the index of the Hamiltonian is well-defined. 

In addition to the chiral invariance we can define particle-hole symmetry. Let us introduce 
the charge conjugation matrix C by 

C(7^)*C-^ = 7^ CiV^yC"^ = -T". (2.3) 

It is always possible to find C for a given set of 7^. Then Hq satisfies 

CHoC-' = -Ho, (2.4) 

where C = CK denotes complex conjugation K followed by the multiplication by C . To 
each particle state if) of an energy £" > we can define a state ipc with energy —E by 

%l)^{x) = C^pix) = Cip*{x). (2.5) 
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Eq. (12.41) ensures the particle-hole symmetry of the spectrum of Hq. 

Before going into detailed analysis of the index of the Hamiltonian we introduce spin((i) 
gauge field 

Mx) = (2.6) 

where 11°'^ = [F", r^]/4 are the spin((i) generators and Aabj are the components of the gauge 

field satisfying A*^^- = Aabj = —Abaj- The gauge symmetric generalization of Hq is then 
given by 

H = -i-fWj + (P{x), (2.7) 
where Dj = dj + Aj is the covariant derivative. For the Higgs field it is given 

Dj4> = V'D^^a, D,4>a = d,4>a + A^^K (2.8) 
The field strength is also defined as usual 

Fij ~S Fabiji Fabij diAabj djAabi ~\~ AaciA(^f)j AacjA^bi- (2-9) 

Since the spin((i) generators commute with 72^+1 and C, the chiral and particle-hole symme- 
tries remain intact by the generalization. 

Unlike Nielsen-Olsen vortex and 't Hooft-Polyakov monopole, which will be treated in 
Sect. [6l we consider the Yang-Mills and Higgs fields as independent classical background. For 
the present we only assume that the gauge potential approaches to pure gauge and the Higgs 
field satisfies 0^ — )■ |0oP at the spatial infinities, where |0oP is a nonvanishing constant. This 
gives rise to a finite energy gap. Due to the particle-hole symmetry mentioned above, any 
eigenstate of nonvanishing energy E is necessarily paired with an eigenstate of energy —E. 
These states are also related with the chiral transformation. The zero modes of if, however, 
can be unpaired. They can be made self-conjugate under the particle-hole symmetry and 
can be regarded as Majorana states. For a general Yang-Mills-Higgs background it is not 
possible to find an explicit form of zero mode wave function. In the chirally symmetric case, 
however, we can investigate the existence or nonexistence of zero modes by computing the 
index of H defined by 

indexif = n+ — n_, (2.10) 

where n± are the numbers of positive and negative 72^+1 chirality zero energy states. The 
index of H is known to be a topological invariant, i.e., it is invariant under continuous 
deformations of the gauge and Higgs fields. The index theorem relates the index with a 
topological invariant of these fields P [101 ITT] . 
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To establish the index theorem we rewrite Eg. (12.101) as 

2 

171 

indexH = lim Tr72rf+i— ^ 

m^O //^ + 

2 



d'^x limtr72,+i ^ 5\x - y), (2.11) 

where Tr stands for the integration over the spatial coordinates as well as the trace on the 
7 matrices. As is well-known, index of the Dirac operator is related with the chiral anomaly 
of the axial current divergence. To show this we introduce an axial current J*(x) by 



J\x) = hm lim tr72.+iy -—— - ^77"^ ' v), (2-12) 

m^o y^x \iH + m iH + M I 

where M is a Pauli-Villars mass to regularize the current at short distances. It is straight- 
forward to show that the axial current divergence can be cast into the form 

d,J\x) = -2 hm lim tr72.+i _^ , - JT^-J7^ ' v)- (2-13) 

At this stage we can take the two limits, m — and M — oo, separately. Eg. (12. lip then 
can be written as 

indexH = -- I dSiJ\x) + Cd, (2.14) 
2 isi-i 

where Cd is the topological index defined by 

and 5*^^ denotes the infinities of the d dimensional euclidean space. In the case of index 
theorem on compact manifolds without boundaries the contribution from the chiral current 
on the rhs of (I2.14p vanishes and the index coincides with q. Since we are working with the 
Hamiltonian (12. 7p defined on a euclidean space, the surface term can be nontrivial [9|[Tm[TT|. 
We expect a nonvanishing contribution to indexif if J^{x) is of order Ixj^''^^ for \x\ — >■ oo. 
To compute the integral on the rhs of Eg. (I2.14p we only need the leading behavior of the 
chiral current at spatial infinities. This will be done in Sect. |H 



3 Topological index 

We have shown that the index can be written as a sum of the topological index (I2.15P and 
the surface integral of the chiral current. The evaluation of the functional trace on the rhs 
of (I2.15P is similar to that of chiral anomalies. In this section we compute q. 
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Eq. f l2.15p can be explicitly written as 



Cd 



M->oo 



lim d X lim 1 172^+1 



M2 



iJ2 + M2 



lim / d'^x 

M-5>oo 



d'^k 



-ikx 



ikx 



Using 



-ikx ^ikx 



—ikx Tj2 Akx 

e ri e 



-iY{Dj + ikj) + (p, 



(3.1) 



(3.2) 
(3.3) 



we can compute the rhs of Eq. (13. ip as 



Cd 



V f ,d [ d^k 
lim / d X / - — -7 

M^oo J J (27l)'^ 



M2 



lim I d'^x 



d'^k ^ 



1 



(27r 



n=0 



72d+l- 



-(A + «A;i)2 + 02 + 



(3.4) 



Due to the presence of 72^+1 the terms with n < d/2 vanish under the trace, whereas the 
terms with n > d/2 do not contribute to the sum in the limit — > 00 as can be easily 
seen by the scaling argument of the momentum variables k — ?■ Mk. This immediately gives 
Cd = in odd dimensions. 

For d = 2N even we obtain 



Cd 



d^x 



-1 



d^^'k 1 /I , . 

(2^ (pttf^'^^^^^^ 



d X6 tr^^jjjj ■ ■ ■ Fi^_^i^^ 



N 



(3.5) 



(27r)^iV! 

where e*^ "*'' is the Levi-Civita symbol in d dimensions. We have also introduced tr^ by 

tr,(---) = tr(r/---), (3.6) 

where 77 is defined by 

(_l)d{d-l)/2 

^ ' -ri---r^. (3.7) 



2'^ 



The overall normalization is chosen to give 



tr^r'^i ■ ■ ■ r"'' = e"^-"". 



(3.8) 
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It satisfies anti-cyclic property 

tr^rr' = -tr^r'r for {rj, T} = 0. (3.9) 

Later we need to compute traces involving 0. In even dimensions 7] anti-commutes with (p, 
whereas it commutes in odd dimensions. 

We have seen that only the gauge field contributes to the topological index and the Higgs 
field is irrelevant to the computation of q. The result (13. 5p coincides with computation of 
index of a Dirac operator without Higgs fields. 



4 Computation of chiral current 

The chiral current defined by (12.121) is written by the regularized fermion propagator. It is 
in general a nonlocal quantity of the background fields. To find its contribution to the index 
of the BdGH we only need the asymptotic behaviors at spatial infinities, where the chiral 
current turns out to become local. In this section we evaluate the asymptotic form of the 
chiral current. 

We first rewrite the current (I2.12p as 



■n^) = / (i^e-tr..,„yi/ - e-'^ (4^1) 



The computation of J*(a;) is similar to the one presented in Sect. |3l We obtain 




|0oP +^2)"+! 

^i^l2d+il\l'k, - i^W, + 0) (^Z7'^fc0 + ll'l'Fki + , (4.2) 

where A is given by 

A = 22A;,A + A'-I0l'+I0or. (4.3) 

At the spatial infinities approaches to a constant |0oP- We further assume 

0^ - l^oT, dj(j), Aj ~ 0(|x|"^) for |x| -> oo. (4.4) 

To find indexif it is only necessary to know the leading 0(|x|~'^"''^) terms of the current for 
|x| — 7- oo. It is easy to convince oneself that the terms of the rhs of Eq. (14. 2 p vanish for 
n < {d—l)/2 because of the trace with 72n+i; whereas the terms with n > d can be ignored 
since they decay faster than |x| for \x\ — )■ oo. Eq. (14.21) then can be written as 



n>{d~l)/2 



xtr72d+i7V(n'A0 + ^7Vi^..-) +0(|a;r'^). (4.5) 
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The k integral can be done for n > (d — l)/2 as 

d'^k 1 T{n+l-d/2) 1 



+1 



(471)^^/%! /i2(n+l)-d- 



(4.6) 



Now the hmits m — and M — )■ oo can be taken safely. Keeping only the nonvanishing 
contributions, we obtain 



T{d/2-l) 



1=0 



{4nY/^{d - I - ly. 



tr72d+i7VSymm 



0(|a; 



where = 0/|0o| and Symm denotes symmetrized product defined by 



SymmA-E- = ^--^ exp[.A + tB] 

n\m\ os^ot"^ 



s=t=0 



Since 'j^Djcj) and 'j^'j'^Fjk are effectively commutative in the trace of Eq. (14. 7p . we can 
simplify the symmetrized product further. We thus arrive at the asymptotic form of the 
chiral current 

P(X) - (-^)'^~^ n2-W I Ip p. . 

(4.9) 

where nonleading contributions are suppressed. The asymptotic chiral current is local with 
respect to the background fields. Unlike the topological index Cd it exists in odd as well as 
even dimensions. 

It is easy to obtain explicit forms of J*(x) in low dimensions. In two dimensions the 
chiral current and topological index are given by 

J\x) = -e'^e'''4>ad,^,-^e'^e'''A,,,, (4.10) 

C2 = d^xe'^e'^'d.Aabj. (4.11) 

These lead to the index 

indexif = — [ dSie'^e^'^^adik- (4-12) 

We see that the topological index is canceled by the gauge field dependent term of the chiral 
current 1101. 
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Similar thing also happens in three dimensions. The topological index vanish identically 
in odd dimensions and the gauge field dependent terms of the chiral current can be cast into 
a total derivative term as 

= -^e'^''e'^Hadikdkk-^e'=h'''>'d^{Aakkk)- (4-13) 
47r 47r 

We thus obtain the index 

indexi/ = ^ I dS,e'^'e''''^adj4>bdkk- (4-14) 

Again the gauge field dependence disappears and the index is only determined by the asymp- 
totic behavior of the Higgs field at the infinities [9]. In the next section we show that this 
holds true in arbitrary dimensions. 



5 Differential geometric approach 

To establish the gauge field independence of indexif it is convenient to introduce Lie algebra 
valued differential forms [16], the gauge potential 1-form and the field strength 2-form, as 



A = ir'^rMxM,,,, (5.1) 

F = ^dx'dx^ Fij = ^dx'dx^ TT^Fabij = dA + A\ (5.2) 
Note that F satisfies Bianchi identity 

DF = dF + AF - FA = 0. (5.3) 
The exterior covariant derivative of (p is given by 

D^ = d4)+[A,$\. (5.4) 

Taking D once more we obtain 

D^^ = [fJ]. (5.5) 

Since the generators of spin((i) commutes with rj defined by fl3.7p . so do A and F in arbitrary 
dimensions. 

The topological index (13. 5p in c? = 2N dimensions can be written in terms of F as 
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where the integral is taken over the entire d dimensional space. If we introduce Chern-Simons 
form ijJ^_i by 

tr,F^ = dul„ (5.7) 
the topological index fl5.6p can be converted to the surface integral 

■u" r ,„ 



where 5*^^ denotes the {d — l)-sphere at the spatial infinities. The Chern-Simons form uj^_i 
can be written as 

= N [ dttTr,A{tdA + t^AY~^- (5.9) 
Jo 

For (i = 2, 4, 6 the Chern-Simons forms are explicitly given by 

ul = tr^ (^A{dAY + ^AMA + ^A^^ . (5.10) 

We now turn to the chiral current (14.91) . The first integral in the rhs of Eq. (I2.14p can 
be written as 



dS,r = / V, (5.11) 

Sd— 1 / crd— 1 

where * J is the Hodge dual of the chiral current 1-form J = Jj(x)da;* and is given by 

*J = jj^e,,,,...,,f'{x)dx''---dx'''. (5.12) 
For the current Eq. (14. 9p . the dual can be written as 

m 

V = 5^ Crf,fctr,0(Z}0)'^-2fc-i^A:^ ^5_^3) 

k=0 

where Cd,k is defined by 
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Taking the exterior derivative of the current and using Eq. fl5.3p and (15.51) . we get 



k=0 

= C,,ktT,{D4>Y-''F' -Y^{d-2k- l)C,,fctr,[F, mm'^^'^'F^' 



k=0 



k=0 



(5.15) 



In deriving this use has been made of the fact that [F, </>] effectively anti-commutes with (j) 
and D(j) in the trace. We can simphfy the trace in the second summand on the rhs of Eq. 
(15151) as 



tr,[F,0]0(D0) 



1 



ii^(t){D(t)Y-^^-^[F''+\. 



k + l 



d-2k~2pk+l 



Inserting this into Eq. (I5.15p . we obtain 

2{d-2k + l] 



d*J 



Cd,k-itT^WY~^^F'' 



k 



(5.16) 



(5.17) 



where use has been made of tr^(D</))°' = and the relation 



_ 2{d-2k + l) ^ 

'^d,k — ^ '^d,k-l, 



k — 1,2, ■ ■ ■ , 



d- 1 



(5.18) 



The former can be verified by noting (f)"^ = 1. 

In odd dimensions the coefficient on the rhs of fl5.17p vanishes. We therefore obtain 



dV = 0. 



This implies that the current can be written as 



2 2'^7r'^(^)! 



(5.19) 



(5.20) 



where Qd-2 is a (c/ — 2)-form. 

In even {d = 2N) dimensions Eq. (I5.17P can be written as 
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(5.21) 



This is reminiscent of chiral anomaly. It also implies that * J can be written as 

-\*J = {-lf^%tT,mr~' + dn,^S,A)-^^y,_,. (5.22) 

Interestingly enough, there appears the the Chern-Simons form. 

It is straightforward to check Eq. fl5.20p and f l5.22p in low dimensions. We give explicit 
expressions for = 2, 3, 4, 5: 

d = 2 : * J = — tr^0d0 tr^A, 

TT TT 

d = ?>: * J = -— tr„0(d0)2 - -dtr„0A, 
Arc Tx 

d = 4: *J = _-Ltr^0(d0)=^ + -^^dtr^ {^^d^A + ]^^A^A^ + -J^tr^ (^AdA + '^A'^ , 
d = b: 

*J = T77r^tr.0(d0)^ + -\dtrJ 0(d0)2A + ^d^A^A + \{^Af + ^AdA + dAA + 
3z7r^ 47r^ I 3 

(5.23) 

We now turn to indexif. We see from Eq. fl5.20p or (15.221) that ^^--2 does not contribute 
to the index and the topological index in even dimensions is canceled by the Chern-Simons 
form. We thus obtain 

indexi/ = (-!)['*'] ^^^^^1^ Ij^mY^'- (5.24) 

The gauge fields disappear completely and the index is determined only by the behaviors 
of the Higgs fields at the spatial infinities. This generalizes the observation of Ref. [10] for 
the two dimensional model of Jackiw and Rossi [13]. In the next section we will introduce 
Yukawa coupling constant g by substituting (f) by gcf) in Eq. (12. 7p . In even dimensions the 
index (I5.24p is not affected by this, whereas an extra overall factor sgii{g) appears on the 
rhs of Eq. (I5.24p in odd dimensions. 



6 Topological configurations in two and three dimen- 
sions 

So far we have considered the Higgs and gauge fields as independent background fields. In 
some model field theories topological objects can be realized dynamically as a solution to the 
field equations. More specifically we assume the gauge-Higgs system in d spatial dimensions 
with a static energy 

^ = /rf'^(^(i^a..,f + ^(A0a)^ + y(|0Or-|0P)^), (6.1) 
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where e is the gauge couphng constant. It becomes stationary for the fields satisfying 

AA0a + A(l - = 0, DiFabij + K^a'^jk = 0, (6.2) 

where k = e^|(/)op and A = Ao|</>oP/2. These nonhnear field equations give rise to topological 
objects, Nielsen- Olesen vortex in two dimensions [1] and 't Hooft-Polyakov monopole in three 
dimensions [2]. We consider the BdG equation with these topological configurations and see 
more closely how the index relation is fulfilled. 

6.1 Vortex in Maxwell-Higgs system 

In two dimensions Eq. fl6.2l) can describe vortices. Ansatz for a vortex with unit vorticity is 
given by 

4)a{x) = h{r)^, Aabi{x) = -eabeij{l - k{r))^, (6.3) 
where h{r) and k{r) are assumed to satisfy the following boundary conditions 

h{Q) = 0, k{0) = l, h{oo) = l, k{oo)=0. (6.4) 
Eqs. (16. 2 p give differential equations for h{r) and k{r) as 

/i" + - = ^ - A(l - h')h, k" -- = Kh'^k. (6.5) 

We see that 1 — h and k decrease exponentially for large r. 

We now turn to indexif for the vortex background. The covariant derivative -Di^a is 
given by 

Di(j)a = h'—^ + eijeabhk—^. (6.6) 

It decays exponentially at r — )■ oo, so does the chiral current *J. We see that indexif 
coincides with the topological index C2 since the chiral current has no contribution to the 
index. It is easy to compute C2- The field strength can be found as 

k' 

Fabij = —^ab^ij — - (6.7) 

Eq. (15.81) then immediately gives indexif = C2 = — 1. 

This might be felt contradictive with the result of Sect. [5] that indexif is determined 
completely by the asymptotic behavior of the Higgs field. It is of course not the case. 
The gauge field is related with the Higgs field by the field equations. In particular the 
chiral current vanishes exponentially as r — )■ oo. Therefore, the contribution from the Higgs 
current, the first term on the rhs of (I4.10p . cancels that from the gauge current, the second 
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term of the same equation, which in turn exactly cancels the topological index of the gauge 
field strength (16. 7p . This implies that the topological invariant (I4.12p coincides with the 
topological index. 

The nonvanishing index suggests that the Hamiltonian given by Eq. (12. 7p has one nega- 
tive chirality zero mode. For the vortex background we can find zero mode of H explicitly. 
To see this we employ the following set of 7 matrices 

7J = a^ ®a\ r^ = (T^®a\ = 1 ® (i = l,2) (6.8) 

The spin (2) generator and chiral matrix is given by 

S3 = = ^ ^ {-if^^^'^T^T'' = 1 ® (6.9) 

Zero mode wave function can be chosen to be chiral. We write it in chiral spinors il)± as 

\v±[x)j 

where x± ^i-re two component eigenspinors with a^x± = =tx±- Each component of the chiral 
zero mode satisfies in polar coordinates 

-ie'^ (dr + -de] ± i^^e'^\u± - g\(t)Q\he^'%± = 



^|0o|/ie^'V + !^-ie^'^ (^dr - ^de^ ± z^^e-'^j v± = 0, (6.11) 

where we have introduced Yukawa coupling constant g by replacing (p with g(j) in Eq. (12. 7p . 
The index (I4.12p is not affected by this change as mentioned in Sect. |5l 

For the negative chirality zero mode we can assume that U- and w_ are independent of 
6. It is easy to check that a normalizable solution is given by 



U- = —isgn{g)v^ = Coexp 







\g4>o\h{r') + LJ^ ) dr' 



(6.12) 



where Co is a normalization constant. In Figure [H we give a plot for the profile of M_(r) 
together with h{r) and k{r). The zero mode wave function is localized around the vortex 
core. 

As for the positive chirality zero mode, we can separate the angle variable by assuming 



where m is an integer, fm and gm satisfy 

d m 1 — k 



dr r 



fm + g\(t>o\hgm = 0, 



d m + 2 1 — k\ ,,,„ 

^ + — —]gm + g\4>o\hfm. = o. (6.14) 
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Figure 1: Profiles of h{r), k{r) and M_(r) witli X = k and g = e/2. 

Tliese lead to the behaviors fm ~ r"^ and Qm ~ r"'""^ as r — > 0. We see that only one of the 
two independent solution is regular at the origin. Such a regular solution, however, contains 
exponentially growing component ~ e^M^ ^ _^ '\Ye thus conclude that there is no 
positive chirality zero mode. 

6.2 't Hooft-Polyakov monopole 

Next we consider a Yang-Mills-Higgs system with spin (3) gauge symmetry in three spatial 
dimensions. The ansatz for the monopole of unit magnetic charge is given by 

U^) = h{r)-, A,U^) = -(1 - A;(r))^^^^^^^^. (6.15) 

Eqs. (16. 2 p are satisfied if h{r) and k{r) obey the following differential equations 

h" + -h' = ^k'^h- X(l-h^)h, k" = Kh^k- ^il-k'^)k. (6.16) 

The boundary conditions for h{r) and k{r) take the same form as Eq. (16.41) for the vortex. 
From Eqs. (16.161) we see the asymptotic behavior k ~ e"^** and 1—h ~ e~^^^ for sufficiently 
large r. No analytic solution is not known for Eqs. (I6.16p . See Ref. [15] for a recent high 
precision numerical study. 

It is straightforward to evaluate indexif. In odd dimensions only chiral current Eq. (14.13^ 
contributes to the index. Note that the covariant derivatives DiCp decays exponentially as 
r — !■ oo and the field strength approaches 

FaM, ^ ^ - j). (6.17) 
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Keeping terms that survive at r — )■ oo, we obtain 

J* ^ -^sgn{g)e'^''e''''4>aF,kbc = -sgn{g)^ (6.18) 
Sir Zirr-^ 

where the overall factor sgn{g) comes from the introduction of the Yukawa coupling constant. 

In odd dimensions the index depends on the sign of g. This immediately gives indexif = 

sgn(5f). As in the vortex case, it is also possible to obtain the same result by computing the 

topological invariant fl4.14p . 

The index obtaind above implies the existence of a zero mode of chirality sgn(5f). For 

the monopole background Eq. (16.151) it is also possible to find the wave function for the zero 

mode. We employ the following representation of the 7 matrices 

Y = a'®l0a\ = l®a''^a'^. (z, a = 1,2, 3) (6.19) 

The spin(3) generators Sq = ^eabc^""^ and the chiral matrix 77 are given by 

Sa = ^1 ® fx'* ® 1, 77 = 1 ® 1 ® (6.20) 



Let us denote the zero mode wave function by chiral components as 
Then ip^ must satisfy 



^ = l"^""!- (6.21) 



(^-ia^ ® Wj + i(T^' ® (T^Aaj ± igl ® (x'^^a^ ip± = 0. (6.22) 

where Aai is defined by Aai = \^abc^bci- These can be cast into 2x2 matrix equations by 
noting {A ® Bilj±)ai3 = Aa'^Bfss{'ip±)-yS = {Aip±B'^)ai3. With this notation Eq. f l6.22p for the 
monopole background Eq. fl6.15p can be expressed as 

-^a^^,<^/±-e,ab{l-k)^a^^±a''T^9\Mhy'^^±^'' = 0. (6.23) 
where are defined by ^± = itp±a^. These have spherically symmetric solutions 

{^±{x))ap = Fi{r)6^p, (6.24) 

where F± satisfy 

F!,{r) = -(^±g\<p,\h{r) + ^-^^^F±{r). (6.25) 

For g > the negative chiral component F_ must vanish, otherwise it grows exponentially 
as r — 7- 00. We thus arrive at the normalizable positive chiral zero mode 

1 , (6.26) 
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Figure 2: Profiles of h{r), k{r) and F^{r) for the monopole background with A = n/2 and 
9 = e/2. 



where is given by 



Co exp 



^7|0o|/^(r') 



A;(r') 



c/r' 



(6.27) 



Again Co is a normalization constant. The zero mode is localized around the monopole and 
the wave function decays exponentially for r — oo. In Figure [2] we give a plot of 
together with h{r) and k{r). The case of < can be analyzed similarly. We obtain one 
normalizable zero mode with negative chirality. This is consistent with the index theorem. 



7 Summary and Discussion 

We have evaluated the index of BdGH of a gauged topological insulator or Yang-Mills-Higgs 
system in arbitrary dimensions by regarding the Higgs and Yang-Mills fields as external 
backgrounds, which can be set up independently. The index can be expressed as a surface 
integral of a gauge invariant chiral current plus topological index of the Yang-Mills fields. 
In odd dimensions the topological index vanishes identically and the gauge field dependent 
terms of the chiral current can be gathered into a total derivative at spatial infinities, giving 
no contribution to the index. In even dimensions the gauge field dependent terms of the 
chiral current can be converted into a total derivative plus the Chern-Simons form, which 
exactly cancels the topological index. We have thus shown that the index of the BdGH is 
determined solely by the asymptotic behavior of the Higgs fields whatever topological charge 
the Yang-Mills field carries. 
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If the behavior of the Yang-Mills and Higgs fields are governed by some effective Hamil- 
tonian, Di(pa and Fij must decay faster than for |x| — )■ oo to ensure the finiteness of 

the Hamiltonian. In such systems nonvanishing topological invariant can be obtained only in 
spatial dimensions less than four. In two dimensions the index of the BdGH is saturated by 
the topological index. In three dimensions only the (pF term of the chiral current contributes 
to the index. This apparently contradicts to the general conclusion that the Yang-Mills field 
does not contribute to the index. It is , however, possible to have expressions for the in- 
dex only in terms of (f) by noting that the gauge fields are related with the Higgs fields by 
Dj4>a = at the spatial infinities. 

We have considered the case of the BdGH containing d order parameters from the be- 
ginning. It is possible to consider other systems with less or more order parameters. The 
evaluation of index of the corresponding BdGH is straightforward. It is rather obvious from 
our explicit calculations that one cannot obtain nontrivial index unless the number of the 
order parameters coincides with the spatial dimensions. Our result is consistent with the 
the topological classification by the Chern number computed from the Berry connection of 
the Bloch wave functions. 
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